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By using a heavy chiral unitary approach, we study the S wave interactions be- 
tween heavy vector meson and light pseudoscalar meson. By searching for poles of the 
unitary scattering amplitudes in the appropriate Riemann sheets , several 1 + heavy 



states are found. In particular, a D*K bound state with a mass of 2.462 ±0.010 GeV 



which should be associated with the recently observed -D s i(2460) state is obtained. 



Q_i! In the same way, a B*K bound state (B s i) with mass of 5.778 ± 0.007 GeV in the 

' bottom sector is predicted. The spectra of the dynamically generated D\ and B\ 



states in the 1 = 1/2 channel are also calculated. One broad state and one narrow 
state are found in both the charmed and bottom sectors. The coupling constants 
and decay widths of the predicted states are further investigated. 
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I. INTRODUCTION 



In recent years, the discovery of various new hadron states stimulated much theoretical 
effort on the hadron spectrum. Among these states, D* (2317) and _D s i(2460) [l] are the 
most attractive ones, because the measured masses are smaller than those predicted in 
terms of most phenomeno logical models 2] (refer to the recent review articles [3]). Many 
physicists presumed that these new states are conventional cs mesons {2, 4], and the others 
believed that they might be exotic meson states, such as tetraquark states molecular 
states , or admixture of cs with molecular component or tetraquark component 

[l(3], and etc.. 

In the corresponding non-strange sector, there are two confirmed D\ states. The narrow 
one is named as £>i(2420) with mass of 2423.4 ±3.1 MeV and width of 25 ±6 MeV, and the 
broad one is entitled as L>i(2430) with mass of 2427 ±26 ±25 MeV and width of 384 + if±74 

nn 

MeV [111, 1121] . Theoretically, the wide-width state with mass of 2325 MeV and narrow- width 
state with mass of 2552 MeV were declared in the x _ SU(3) approach [71] . In Ref. 8], these 
two states were considered as quasi-bound states and were used to determine the unknown 
coupling constants in the next-to-leading order heavy chiral Lagrangian. The reproduced 
masses and widths are M Dl = 2422 MeV and T Dl = 23 MeV and _M p > = 2300 MeV and 

= 300 MeV, respectively. In many other references [13I . . 
were proposed as the conventional cq states, for instance, in Ref. 
as the mixed 1 P\ and 3 Pi cq states with a mixing angle of 
measured widths. 

Since the _D* (2317) and _D* 1 (2460) have the same quantum number J, the isospin, and 
S, the strangeness, except J, the total angular momentum, in the S wave DK and D*K 
channels, respectively, it would be constructive to study the DK and D*K interactions so 
that one can see whether these two states could have similar molecular structure. In the light 
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171 ]. these two states 



151 ]. they were deliberated 
a 35° obtained by fitting 



hadron system, the chiral unitary approach (ChUAJhas achieved 
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great success in explaining 
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23j, 
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251 ] . Some well 



the meson-meson and meson-baryon interactions 
known hadrons can be dynamically generated as the quasi-bound states of two mesons or 
a meson and a baryon 



m 



20. 
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24j . for instance, the lowest scalar states a, / (980), a (980), k 



25], and etc.. Then, the approach was extended to the heavy hadron system, 



called heavy chiral unitary approach (HChUA) 26]. In terms of such an approach, the S 
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wave interaction between the heavy meson and light pseudoscalar meson was studied, and 
some bound states and resonances were predicted, for example, the _D* (2317) state as a 
DK bound state at 2.312 ± 0.041 GeV, a BK bound state at 5.725 ± 0.039 GeV |2Q|. In 
the same approach, D$ and Bq in the (J, S) = (1/2, 0) channel were also investigated. 
As a result, one broad state and one narrow state were predicted in both the charmed and 
bottom sectors [26 ]. 

In this paper, we study the S wave interaction between heavy vector meson and light 
pseudoscalar meson and search for J p = 1 + heavy mesons in both strange and non- strange 
sectors. The couplings of various coupled channels to the generated states and the decay 
width of the isospin symmetry violating process D s i(2460) + — > D* + n° are studied as well. 



II. DYNAMICALLY GENERATED HEAVY-LIGHT 1+ STATES 



To describe the interaction between the heavy vector meson and the light pseudoscalar 

□ n 

meson, we employ a leading order heavy chiral Lagrangian [3, 127J 

1 



(1) 



where f n = 92.4 MeV is the pion decay constant, P* represents the heavy vector mesons 
with quark contents (Qu, Qd, Qs), namely (D*°, D*+, D* + ) and (B*~, B*°, B*) for the 
charmed and bottom sectors, respectively, and $ denotes the octet Goldstone bosons in the 
3x3 matrix form 



7Y~ 



7T 



K+ \ 



V 



(2) 



We are interested in (J, S) = (0, 1) and (J, S) = (1/2, 0) systems. Usually, such a 
system can be characterized by its own isospin. Based on the Lagrangian in Eq. 
amplitude with a definite isospin can be written as 



the 



V i I j {s,t,u)e ■ e' 



r 1 



U £ • £ 



(3) 



where i and j denote the initial state and the final state, respectively, s, t, u are usual 
Mandelstam variables, and e and e' are polarization vectors of the vector states in the initial 
and final states, respectively. In the 1 = case, there are two coupled channels. The channel 
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label 2 = 1 and 2 specify the D*K (B*K) and D*r) [B*r]) channels in the charmed (bottom) 
sector, respectively. In the 1=1/2 case, three coupled channels exist. The channel label 
% = 1, 2 and 3 in this case denote the D*n (B*ir), D*r] (B*rj) and D*K (B*K) in the 
charmed (bottom) sector, respectively. The corresponding coefficients are tabulated in 
Table [H It should be mentioned that, in the coupled channel calculation, the thresholds 
of the channels with the light vector meson and the heavy pseudoscalar one are relatively 
higher than those with the light pseudoscalar meson and the heavy vector one. For instance, 
in the charmed 1 = 1/2 case, the lightest combination of a light vector meson and heavy 
pseudoscalar meson is p + D, and the sum of their masses are almost forty MeV heavier than 
that of the heaviest combination of a light pseudoscalar meson and heavy vector meson, say 
K + D*. Thus, in the concerned energy region near the thresholds of the channels with later 
combinations, the contributions from the channels with former combinations are expected 
to be less important, and hence can be neglected for simplicity. 



TABLE I: Coefficients CL in Eq. ©. 



r<o rfQ r<o 

u ll u 12 u 22 


n l/2 r l/2 r l/2 r l/2 
u ll u 12 u 22 u 13 


^1/2 
°23 


^1/2 
°33 


-2^0 


-2 


2 


-1 



In ChUA, the unitary scattering amplitude can be expressed by the algebraic Bethe- 
Salpeter equation [3]. The full unitary amplitude for the S wave scattering of vector and 
light pseudoscalar mesons can be written as [3, Q] 

T T (s) = -[1 + V I (s)G(s)(l + ^)]-V^), (4) 

where the polarization vectors are dropped because they are irrelevant to the pole searching. 

In the equation, My is the mass of the vector meson in the meson loop and q on represents the 

on-shell three-momentum in the center of mass frame. V^s) is in the matrix form with its 

elements being the S wave projections of V/j(s,t,u). The non-zero element of the diagonal 

matrix G(s) is the two- meson loop integral 

f d A q 1 1 

Gll{s) = 1 J (2vr) 4 g 2 - ml + is {p 1 + p 2 - qf — My + ie ' (5) 

where is the mass of the light pseudoscalar meson in the loop. The loop integral is 
calculated in terms of the dispersion relation with a pre-selected subtraction constant a(/x) 
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2l| . The subtraction constant can be fixed by matching the calculated loop integral with 



that calculated by using the three- momentum cut-off method [26J. The matching point is 
taken at Md* + m K for the charmed sector and M B * + m K for the bottom sector, respec- 
tively, because we are interested in the energy region around the threshold. With the same 



consideration shown in Ref . 26] , the estimated three-momentum cutoff q ma x is in the region 
of 0.8 ± 0.2 GeV. The corresponding values of a(/x) and q m ax are tabulated in Table HT1 The 
resultant loop integration curves in two different methods are plotted in Fig. [TJ It is shown 
that they are very similar with each other in the region around the matching point. With 
the estimated a(/z) values, the unitary scattering amplitudes can be calculated. 

TABLE II: Corresponding a(/i) values and q ma x values with /j, = mp for the charmed sector and 
H = tub for the bottom sector, respectively. 



Qmax (GeV) 


0.6 


0.8 


1.0 


a(m D ) 


-0.639 


-0.714 


-0.752 


a(m B ) 


-0.0764 


-0.101 


-0.113 
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FIG. 1: The real parts of the loop integrals calculated by using the cut-off method (dashed lines) 
and the dispersion relation method (solid lines), (a) D*K loop, (b) B*K loop. 



The poles of the full scattering amplitudes in the 1 = 0, S = 1 channel in both the 
charmed sector and bottom sector are searched first. It is shown that on the first Riemann 
sheet of the energy plane, there is a pole located in the real axis below the lowest threshold of 
the coupled channels in either the charmed sector or the bottom sector. The resultant pole 
positions with different a(/z), which correspond to q max = 0.6, 0.8 and 1.0 GeV , are tabulated 
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TABLE III: Poles in the (I, S) = (0, 1) case. 



Qmax 


(GeV) 


0.6 


0.8 


1.0 




(GeV) 


2.472 


2.459 


2.452 


B s i 


(GeV) 


5.785 


5.775 


5.771 



in Table respectively. These poles are apparently associated with the D*K bound 
state and the B*K bound state, respectively. More specifically, when a(rno) = —0.714, 
corresponding to q max = 0.8 GeV, the mass of the D*K state (Z? s i) is about 2459 MeV, 



which is quite consistent with the measured value of the _D s i(2460) state [11] 

M Dsli24m) = 2458.9 ± 0.9 MeV. (6) 

Taking into account the uncertainty of the subtraction constant, the resultant masses of the 
D s i state and the undiscovered B*K bound state, namely B s i, are 

M Dsl = 2.462 ±0.010 GeV, 

M Bal = 5.778 ± 0.007 GeV. (7) 

The predicted mass of the B*K bound state is consistent with the simple estimate of 5778 
MeV in Ref. [9]. 

In order to show the effect of the coupled D s r] (B s rj) channel explicitly, we also present 
the single D*K (B*K) channel result in Table HVl It is shown that the single channel results 
are slightly larger than the coupled channel results. It implies that D sl (2460) (B s i) can be 
regarded as a bound state of D*K (B*K) with a tiny component of D s r\ (B s rj). 

TABLE IV: Poles in the (J, S) = (0, 1) case in the single channel approximation. 



Qmax 


(GeV) 


0.6 


0.8 


1.0 




(GeV) 


2.478 


2.467 


2.462 


B s i 


(GeV) 


5.786 


5.779 


5.775 



In the / = 1/2, 5 = case, the poles are located on nonphysical Riemann sheets. Usually, 
for a certain energy if Imp cm is negative in all the opened channels, the pole obtained would 
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correspond more closely with the physical one. There are two poles in the charmed (bottom) 
sector. The width of the lower pole is broad and the width of the higher one is relatively 
narrow. We tabulated these results in Table [V] In the charmed (bottom) sector, the lower 





TABLE V: Poles in 


the (I, S) 


= (i 0) 


case. 




Qrnaz 


. (GeV) 


0.6 


Oi 


5 


1.0 




Di 


(GeV) 


2.245 -iO. 106 
2.599 - i0.043 


2.239 - 
2.585 - 


i0.094 
i0.044 


2.236 - 
2.578 - 


i0.088 
i0.043 


Bi 


(GeV) 


5.586 - i0.118 
5.884 - i0.027 


5.579 - 
5.875 - 


i0.108 
i0.025 


5.576 - 
5.870 - 


i0.103 
i0.023 



pole is located on the second Riemann sheet (Imp cm i < 0, Imp cm2 > 0, Imp cm3 > 0, where 
p cm i is the momentum of one of the interacting mesons in the center of mass frame in the i-th 
channel) and should be associated with the D*ir {B*ir) resonance. This state should easily 
decay into D*ir (B*tt). The higher pole in the charmed (bottom) sector is found on the 
third Riemann sheet (Imp cml < 0, lmp cm2 < 0, lmp cm3 > 0) and should be associated with 
an "unstable" D*K (B*K) bound state due to its narrow width. It should be mentioned 
that the pole structures of 1 + states here are very similar to those of + states 26J], but are 
different from that of the /o(980) state where only one pole located on the second Riemann 
sheet and one shadow pole on the third Riemann sheet |28|] . The origin of the difference 
comes from the fact that there are two coupled channels, i.e. the nir and KK channels, in 
the /o(980) state case, while there are three coupled channels in the 1 + state case. 

The fact that two poles in the different Riemann sheets should be associated with two 
different 1 + states in the I = 1/2 channel can be confirmed by checking the curve structure of 
the absolute values of the unitary scattering amplitudes in the coupled channel case, because 
such curve structure is closely related to the structure in the corresponding invariant mass 
spectrum. The absolute values of the unitary scattering amplitudes in the coupled-channel 
case for the D*n — > D*tc and D*K — > D*tt processes in the I = 1/2 channel are plotted in 
Figs. [2] (a) and (b), respectively. It is shown that there are one broad peak and one narrow 
dip in the D*n — > D*tt scattering amplitude and two peaks in the D*K — > D*n case. In 
the bottom sector, the curve structure of scattering amplitudes is the same and will not be 
demonstrated here for simplicity. 
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FIG. 2: The absolute values of unitary scattering amplitudes in the coupled-channel case for the 
D*tt -» L>*vr and Z)*^ -> D*tt processes in the 1 = 1/2 channel, (a) D*tt -» L>*vr, (b) L>*i^ -» L>*vr. 



Considering the uncertainty of a(/i) in Table [VT we predict the masses and the widths of 
the broad Di and Bi states as 

M Dl = 2.240 ± 0.005 GeV, T Dl = 0.194 ± 0.019 GeV, 

M Bl = 5.581 ± 0.005 GeV, r Bl = 0.220 ± 0.015 GeV, (8) 

respectively, and the masses and the widths of the narrow D\ and B\ states as 

M Dl = 2.588 ± 0.010 GeV, T Dl = 0.087 ± 0.001 GeV, 

M Bl = 5.877 ± 0.007 GeV, r Bl = 0.050 ± 0.003 GeV, (9) 

respectively. 

These results mean that if we believe that the -D* (2317) and _D* 1 (2460) states in the 
(/, S) — (0, 1) channel are really dominated by the molecular state structure, the predicted 
two 1 + states in the (I, S) = (1/2, 0) channel should also exist. As mentioned in the In- 
troduction, the experimentally established 1 + states _Di(2420) and _Di(2430) are compatible 
with the conventional cq interpretation. It implies that there are no experimentally estab- 
lished states as the candidates for the predicted (/, S) = (1/2, 0) states. Why these two 
states have not been observed are complex. At the present time, it is not at all clear about 
what to expect with regards to production of these states. In fact, finding a new state does 
not only depend on its high production rate, but also relate to, in a large extent, the data 
measurement and analysis, which are usually affected by many factors, for instance the data 
statistics, the background, the width of the state, the complexity of the spectrum structure 
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in the vicinity of the state, and the suitable channels for producing and detecting such a 
state, and etc.. For instance, one of the possible reasons which makes the observation of 
the predicted states difficult is that the production of such states would be suppressed with 
respect to conventional ones. This is because: (1) The predicted two D\ states are quasi- 
bound states of other two mesons. From the viewpoint of quark degrees of freedom, there 
should be at least four quarks in the Fock space. Therefore, the production of such states 
would be suppressed in comparison with producing a cq state due to the necessary creation 
of an additional quark anti-quark pair. (2) The widths of predicted D\ states are comparable 
with those of the corresponding conventional cq states, so that the couplings between these 
states to the D*tt state would be similar with the cq states. In other words, the signals of 
the predicted D\ states would be suppressed compared with those of the conventional cq D\ 
states in the D*tt spectrum, and a definite observation of such states becomes difficult. It 
seems that an even larger data set with higher statistics and further careful data analysis, 
as well as theoretical study are necessary. 

Moreover, it is interesting to compare our predicted states with the quark model predicted 



conventional qq states shown in Table IVII 13L llo . l29j ]. We find that the lowest masses of 
the conventional Qs states with J p = 1 + are larger than those of our quasibound D s i and 
B s i states, respectively, and the lowest masses of the conventional Qq (q = u, d) states with 
J p = 1 + are sited between the masses of our lower and higher predicted states. Hence, if 
one can find a state with a mass much lower than the quark model prediction, it might be 
the lower state in our prediction. 



III. COUPLING CONSTANTS AND DECAY WIDTHS 

The nature of dynamically generated states can be further studied by calculating the 
coupling constants between these states and the particles in the coupled channels. The 
coupling constants are related to the Laurent expansions of unitary scattering amplitudes 
around the pole [30J 

Tij = -Jtil \- 7o + 7l ( s - s po i e ) H , (10) 

■S Sp i e 

where gi and gj are the coupling constants of the generated state to the i-th and j-th chan- 
nels. The product gigj can be obtained by calculating the residue of the unitary scattering 
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TABLE VI: Masses of l 1 ^! and l 3 Pi heavy-light mesons in quark models. 



Ref. [13] Ref. [29] Ref. [16] 





2.44 


2.46 


2.490 


Di(l 3 P^ 


2.49 


2.47 


2.417 




2.53 


2.55 


2.605 


D sl (l 3 Px) 


2.57 


2.55 


2.535 


B^Px) 




5.78 


5.742 


B 1 (1 3 P 1 ) 




5.78 


5.700 


B sl (l l Pi) 




5.86 


5.842 


B s i(l 3 Pi) 




5.86 


5.805 



amplitude at the pole [21] 



(11) 



giQj = lim (s - s po i e )Tij. 

As a typical example, we calculate the coupling constants for the _D s i(2460) state with 
a subtraction constant which corresponds to q ma x = 0.8 GeV, because under this condition 
the empirical mass of _D s i(2460) can be excellently reproduced. The resultant coupling 
constants are tabulated in Table IVHI for the states with 1 = and Table IVIIII for the states 
with 1=1/2. 

Comparing the coupling constants in Tables IVHI and IVIIII with theose of the generated 
scalar states in Tables VI and VII in Ref. 26| . one finds that the values of corresponding 
coupling constants are close. This manifests the heavy flavor spin symmetry 311] respected 
in the leading order heavy chiral Lagrangian 22J . In Table IVHI the data also show that g\ 



TABLE VII: Coupling constants of the generated D s \ and B s \ states to relevant coupled channels. 
In this case, g\ and g% are real. All units are in GeV. 





Masses 


M 


\92 1 


D sl 


2.459 


10.762 


6.170 


B s \ 


5.775 


23.572 


13.326 
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TABLE VIII: Coupling constants of the generated D\ and B\ states to relevant coupled channels. 
All units are in GeV. 





Poles 


9i 


\9i\ 


92 


\92 1 


93 Iflfcl 




2.239 


- i0.094 


8.157 + i5.135 


9.639 


-0.202 + 


i0.059 0.211 


4.919 + ^3.053 5.790 


Di 


2.585 


- i0.044 


0.145 + i3.306 


3.309 


-6.893 - 


i2.237 7.247 


-11.060 + ^1.165 11.121 


Bi 


5.579 


-*0.108 21.439 + ill.861 24.502 


-2.222 - 


i0.752 2.346 


13.517 + Z6.906 15.179 


Bi 


5.875 


- i0.025 


0.295 + i6.619 


6.626 


-14.553 - 


- i4.892 15.353 


-24.759 - i0.874 24.775 



is larger than g 2 in the charmed (bottom) sector. This indicates that the coupling between 
the D s i (B s i) state and the states in the D*K (B*K) channel is larger than that between 
the D s i (B s i) state and the states in the D*r] (B*rj) channel. It also reflects the fact that 
the generated state is a D*K (B*K) bound state. In Table IVIIIt the largest coupling 
constant \g±\ for the lower state is associated with the D*n (B*tt) channel and the largest 
coupling constant \g^\ for the higher state is connected with the D*K (B*K) channel. This 
is consistent with our finding in the pole analysis, namely the lower state is a D*tt (B*tt) 
resonance and the higher state associates with a D*K (B*K) quasi-bound state . 

The coupling constants also show that the largest component of the lower D\ state is 
D*7r whose quark contents are cnnn, where n denotes the u or d quark, and the largest 
component of the higher Di state is D*K whose quark contents are csss. On the other hand, 
the _D s i(2460) state, in principle, is a D*K bound state, and consequently, the dominant 
quark contents of this state are cnns. Thus, from the quark contents of these states, one 
can expect that the mass of D s \ should be a value between the masses of the two D\ states. 
The same qualitative statement can be applied to the bottom sector. 

Furthermore, the dynamic nature of a state can be characterized by its decay properties. 
Since the quark contents of predicted molecular states are different from those of conventional 
qq states, their decay properties are expected to be different. 

In order to estimate the order of magnitude of the decay width of _D s i(2460), we calculate 
the decay widths of the isospin symmetry violating decay processes -D s i(2460) + — > D* + tt° 
and I?gi(5775) — > B*°7r° can be estimated through ir -r] mixing [32] shown in Fig. [3j By 
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D s1 (2460) + 




FIG. 3: Isospin symmetry violating decay process £> s i(2460) + — > Z)* + 7r° via 7r°-?7 mixing 



using the formula 



Pcm ^^ l 92U v e x {p) ■ e x '*{p LJ 2 



where M is the mass of the initial meson, p cm denotes the three-momentum in the center of 
mass frame, e x (p) and e x ' (pi) represent the polarization vectors of the -D sl (2460) (_B s i(5775)) 
and D* (B*) states, respectively, and J2x Sa' describes the sum over the final states and 
the average over the initial states. By using the 7i°-rj mixing amplitude t nv = —0.003 GeV 
obtained from the Dashen's theorem [33j , we get the decay widths 

r( J D sl (2460) + -> D* s + 7i°) = 11.41 keV, 
r(fi sl (5775)° -»■ B* s \°) = 10.36 keV. (13) 

It should be mentioned that, the explicit isospin breaking term is not the only source of 
isospin breaking, the mass difference of isospin multiplets has already generated such a 
breaking. Thus, the given result can only serve an estimate of the order of magnitudes of 
the widths. 

The decay properties of predicted D\ {B\) state can also be briefly estimated. For the 
higher D\ [B\) states, it can strongly decay into the two opened channels: D*tt and D*rj 
(B*it and B*i]). Although the channel with 77 have much smaller phase space than that with 
7T, due to the relatively larger contribution from hidden strangeness, which can be seen from 
Table EH!! the branching fraction of the D\ [B\) state decaying into the final state with rj is 
comparable with or even larger than that with it. The ratio of these two branching fractions 
can roughly be estimated by using corresponding coupling constants given in Table IVIIII 

r( jDl p*rj) 
T(D 1 -> D*tt) 
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As to the conventional Di (Bi) state which has a mass larger than the D*r] (B*tj) threshold, 
the branching fraction of decaying into D*rj (B*rj) would be much smaller than that into 
D*n (B*ir) due to both OZI suppression [34j and phase space suppression. Thus, it would be 
easy to distinguish this state from the higher D\ (B\) state by measuring the ratio defined 
in Eq. (1111) ■ For the lower D\ (B\) state predicted here, it has only one opened channel 
D*tt (B*7r). So it would be difficult to distinguish the lower state from the conventional Di 
(Bi) state with similar mass and width by using strong decay properties. Yet, the radiative 
decays of such states might be different. The concrete consequences should be investigated 
in future. 



IV. CONCLUSION 



We study the dynamically generated axial heavy mesons which have the same quantum 
numbers of the conventional cs and cq states in the (I, S) = (0, 1), (I, S) = (1/2, 0) systems 
in the framework of coupled-channel HChUA. In the (J, S) = (0, 1) and J p = 1 + system, 
there are two coupled channels: D*K and D*r] for the charmed sector and B*K and B*i] 
for the bottom sector, respectively. In the coupled channel calculation, the channels with 
a light vector meson and a heavy pseudoscalar meson are not considered due to their less 
importance. 

By searching for the pole of the unitary coupled-channel scattering amplitude on appro- 
priate Riemann sheets of the energy plane, we find a state with mass of 2.462 ±0.010 GeV in 
the (J, S) = (0, 1) system. This state should be a D*K bound state with a tiny D*rj compo- 
nent. We would interpret it as the recently observed charmed meson D sl (2460). In the same 
way, we predict a B*K bound state with mass of 5.778 ± 0.007 GeV in the bottom sector. 
In the (I,S) = (1/2,0), J p = 1 + system, there are three coupled channels: D*tt, D*rj and 
D*K in the charmed sector and B*ir, B*r] and B*K in the bottom sector, respectively. We 
find two poles in the nonphysical Riemann sheets in both the charmed and bottom sectors. 
In the charmed (bottom) sector, the lower pole is located at (2.240 ± 0.005 — iO. 097 ± 0.009) 
GeV ((5.581±0.005-i0.110±0.007) GeV). The state associated with this pole will strongly 
decay to D*n (B*tt) and have the largest coupling with the decayed particles. Therefore, 
this pole should be associated with the D*ir (B*7r) resonance. The higher pole is positioned 
at (2.588 ± 0.010 - i0.043 ± 0.001) GeV ((5.877 ± 0.007 - i0.025 ± 0.002) GeV). The state 
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associated with this pole has two decay channels D*n and D*r] (B*tt and B*r]) and has 
the largest coupling with decayed particles in the D*K (B*K) channel. Thus, this pole 
should be associated with a quasi-bound state of D*K (B*K). If one believes that the 
corresponding states in the (I, S) = (0, 1) and (/, S) = (1/2, 0) systems have similar 
S wave molecular state structures, the two predicted states should exist. The estimated 
order of magnitudes for the widths of the _D s i(2460) + and S s i(5775)° states is about 10 
keV. The decay properties of predicted Di (Bi) states are also briefly discussed. Study the 
channels where the final states include D*t] and D*n (B*r] and B*tt) would be helpful to 
find predicted Di (Bi) states. 
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